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ABSTRACT 


By Fourier-transforming the author* s recently proposed state 
functional formalism for the BBGKY hierarchy, a new perspective of 
non-equilibrium statistical mechanics is given: the basic equation 

is formally very close to the Fokker-Planck equation and may readily 
be modified to a universal master equation (with irreversibility) by a 
slight change. Hence the problem reduces to one of a generalized 
random-walk such that the stochastic quantity to be considered is the 
particle-number density in the one-body phase space. A general 
solution is formulated for the weak interaction case. 


1. Introduction 


Recently the functional formalism of classical non-equilibrium 
statistical mechanics was proposed by the author 1 to search for general 
closed-form solutions for distribution functions governed by the BBGKY 
hierarchy for the infinitely-many-particle system with a finite average 
density. There the state functional, which embraces all distribution 
functions in itself as the coefficients of the functional Taylor ex- 
pansion (except for the numerical factors), was introduced. It is 

governed by a single second-order functional-differential equation. This 

2 

formalism originates with Bogoliubov but a basic modification was neces- 
sary for a meaningful functional calculus to be performed within the 
formalism. Be that as it may, it would be quite inconvenient if the 
state functional is too fine-grained to represent a real irreversible 
process, even if the solution could be numerically pursued. In fact, 
on the basis of the exact dynamics, i.e., the Liouville equation or its 
equivalent, we can hardly expect all possible arbitrary states to de- 
velop into a unique, steady (e.g. equilibruim) state. This is obvious, 
if we remember the time-reversibility of the equation as well as the fact 
that the entropy in Gibbs's sense is conserved so that the maximum entropy 
is never achieved, starting from a state of lower entropy. In order to 
establish an entropy-productive non-equilibruim statistical mechanics, 
some modification of the basic equation, comparable with coarse -graining 
procedure, seems to be unavoidable in our state functional formalism also. 
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The first aim of this paper is to find the best irreversible per- 

3 

spective for the evolution of the state functional' „ We first notice that 
the Fourire-transformed basic equation in the cylinder functional approach 
(established in Sec. 2) has a structure very close to the Fokker-Planck 
equation which describes a Markov stochastic process, but retaining the 
time-reversibility as a natural consequence o^ its equivalence to the 
BBGKY hierarchy (in the limit when the number of particles tends to in- 
finity). The best method of modifying the equation into a real Fokker- 
Planck equation at the expense of losing the time-reversibility is 
described in Sec. 3 * As a result, our problem reduces to nothing but the 
problem of a generalized random-walk played within the special function 
space equivalent to a certain Riemannian space by the so-called stochastic 
particle-number density function in the one-body phase space. The 
mathematical procedure of such a modification is unique in the sense that 
any other Markov process expected to simulate non-equilibrium statistical 
mechanics is not closer tc the dynamical process than the present one; 
but it seems somewhat difficult to re-express the procedure in simple 
conventional physical words, such as coarse -graining . 

As the paper ! s second object, the general method is applied to the 
case of a gas or a weakly interacting plasma in Sec. I 4 -. In this case the 
equivalent Riemannian space reduces approximately to a Euclidean space, so 
that the complexity in calculation is reduced. The closed-form general 
solution of the initial -value problem is formulated in terms of a repeated 
multiple integral; which can in practice be conveniently solved by the 
Monte Carlo quadrature. Since in our theory the accuracy of calculation 
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is increased, simply by increasing the multiplicity M of cylinder functional 

but not by changing (or deepening) the formulation in a complicated way 

and, moreover, since there is no difficulty in principle in applying the 

theory to the strong interaction case, the proposed approach seems to have 

the proper advantage which is not found in other analytical approaches to 

the BBGKY hierarchy or the Liouville equation which exploit the more or 

2 If 5 

less sophisticated series expansion in a small parameter, * ’ which may 
not necessarily be convergent. 

Another type of statistical-dynamical research, which may be also cal- 
led a functional approach, has been developed by considering the evolution 
of the probability on an ensemble of the events, a member of which evolves 

strictly according to some dynamics such as the Vlasov equation, the 

6 7 

Langevin equation, etc. ’ Of course, this research is different in 
principle from ours; except for the case -f taking the Klimontovich equation 
as the dynamics, when the characteristic functional equation is formally 
equivalent to the BBGKY hierarchy, as was verified by Nakayama and Dawson . 
However, it will be seen (in Sec. 5) that if the Vlasov equation is taken 
as the dynamics, there is a somewhat close relation between this ensemble 
mechanics and the presented perspective of our formalism. 

The discussion is restricted to a single -component system for sim- 
plicity. Extension to a multi -component system is straightforward^. 
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2. Fourier-transformed State Functional Formalism 

According to the previous paper * if the s-body generic distribution 
function is denoted by F g 5 the state functional is defined as 


. s 


*(y. t) = 1 + l — f • • / f s (x 

S=1 s * X X 


x s , t)y(x 1 )* •y(x s )dx 1 **dx s , (2. l) 


where i = /-l, X i3 the- entire one-bcdv phase space, x n , --x ( X, t is 
the time variable and y(x) is a real-valued function. Once the state func- 
tional is introduced, the BBGKY hierarchy for the inf initely -many -particle 
system with a finite average density is completely replaced by the follow- 
ing single basic equation with functional differentiation: 


= if y(x)[H (x) ; Jdx - i / f (y(x)y(x ' ) - iny(x) - iny(x' )} 

dt ' ^ x i6y(xj c. 

* [♦(k - q' | ); -= ^ ]dxdx' (2.2) 

i 6y(x)6y(x') 


where H^x) is the one-body Hamiltonian, [ ; ] the Poisson bracket, n the 

average number density and 4>(|q - q. ’ | ) the interaction potential between 

particles in which q denotes the displacement vector in the physical space. 

9 10 

The operator 6/6y(x) denotes a functional derivative.’ \p should be sub- 

l 

Ject to the following two associative conditions: 


*(0, t) ■ 1 ; (2. 3) 

lim / — — dx/V = -J/ , (2. 1) 

V-*-» *X i6y(x) 

where V is the volume containing the entire system. Both conditions 

emerge from the form of (2. l) on using the definition of F . 

s 
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I’ 


Here we introduce a kind of cylinder functional approach to func- 
tionals. If {s.(x)} is an orthonormal function set in X, we can calcu- 
J 

late a = / s (x)y(x)dx to define a new function: 

J X J 


y^(x) = l a,s (x) 

J-l J J 


( 2 . 5 ) 


M M 

We call \p(y n ) the cylinder functional \p l (y), which is still a functional 

of y through the definition of a^, though it can essentially be considered 

a function of M variables {a }. Then ip is understood as in the limit 

J 

M + » if it is convergent (for example with respect to the "maximum" 
norm). If > is replaced by ip , the basic equations (2. 2) - (2. U) be- 
come more feasible to treat, since then we have the relation 


x M 6a, a M 

— - 1 


= y s . (x ^ — 

<$y(x) 6yfx) 6a^ j=i J 


( 2 . 6 ) 


and thus, (2. 2) becomes just a partial differential equation. From 

a practical point of view we always consider ^ first, assuming that the 

r M i 

converging sequence } of state cylinder functionals exists. Then wc 

M / \ M 

can introduce a new functional p (x,t; as the Fourier component of tp ; 

namely we have 


V M (y, t) = / 0 M (z, t)exp{i/y‘^(x)z M (x)dx}6z M . 




(2. 7) 


where 


M 


/ 6z M =["••/" n db /(2 it) ' 
A M — J-l J 


1/2 


( 2 . 6 ) 


M r 

andz (x) = 2,hjS^(x). For simplicity we omit all the superscript M 
hereafter, unless necessary for clarity, keeping in mind that every 
functional is a cylinder functional. 


.sttiMkl 
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iM a result of the transformation (2. 7) > the basic equation (2. 2) 


transforms to the following equation for p : 


3p_ _ _ f 

* 


6z(x) 


z(x)p}dx + i / / 


2 'x X | 5z(x)6z(x' ) 


[[♦( k - a' I ) ; 


z(x)z(x' ) ]p}dxdx' , 


and Q denotes the nonlinear operator such that 

Qz(x) = [H (x); z ( x ) ] + n / dx'[<j)(|q - q'|); z(x)z(x')] . (2.10) 

1 X 

It is noted 1 that 3z/3t = Qz is nothing but the self-consistent Vlasov 
equation and therefore Q may be called the Vlasov operator. Corresponding 
to (2. 3) and (2. 4), we have the conditions: 

/ pfiz = 1 ; (2. 11) 

A 

/ p{lim / z(x)dx/V - l}exp{ij y(x)z(x)dx} 6z * 0 . (2. 12) 

A V-*® X X 

From (2. 12) p should vanish for all z(x) unless z(x) is such that 


lim / z(x)dx/V = 1 , (2. 13' 

V-«° X 

M 

which defines a hyperplane in R . We may call this the hyperplane A^. 
The restriction (2. 12) is conservative in time, since the conservation 
of (2. 4) was proved by the previous paper. Therefore, without any in- 
consistency we have only to consider the basic equation (2. 9) always 
in this hyperplane rather than R^. Accordingly A in (2. 7) and (2. li.) 

may be read as A on the understanding that f 6 z indicates the volume 

A 1 

integral over A^. The time-invariance of the condition (2. 11) is evi- 
dent by virtue of (2. 9). 


ililMWiilWl+l: 
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hierarchy or the Liouville equation in the limit N = nV ■+ « (n: finite), 
any dynamical property of the particle system is preserved in our 
formalism. It is easy to verify the time-reversibility of the basic 
equation (2. 9) ; if p(z(q,p), t) is a solution, then p(z(q,-p), -t) be- 
comes another solution, p denoting the momentum part of the vector x. 

Also it is possible to derive the conservation laws of mass, momentum, 
angular momentum and energy directly on the basis of (2. 9)« In this 
case, on comparing (2. 1 ' with the series expression of (2. 7) in powers 
of y, we first have 

F^Cx) = J ^ p(~, t)z(;c)Sz , (2. 14 ) 

F 2 (x, x') = / A p ( z , t)z(x)z(x' )6z (2. 15) 

and then can proceed with examining the time -evolution of each physical 
quantity related with p through the above relations. 

As a representative example let us examine the momentum case. From 
(2. 14) and (2. 9), we calculate 

— / npF (x)dx = n / J pz(x)dx6z 

9t X 1 A 1 9t X 

= n / j pQz(x)pdx6z 

A 1 X 

where use was made of partial integration, keeping in mind that |p| 0 

2 

at infinity (/ z dx * °°) in A-. in order for f A |p|6z to exist. As is easily 
X 1 ' A 1 

seen, there is no effect of the term with a second-order functional-deriva- 
tive in (2. 9) in this case, since the integral over X in the left-hand side 
is not a functional higher-order than linear in z. A further calculation 


gives 
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~ / npF dx = n / p {/ p[H ; z(x)]dx + n/ J p[ <p ; z(x)z(x' ) ]dxdx' }5z 
3t X 1 A 1 X 1 XX 

3H 

= n / p { — / — — z(x)dx - nf f (— + — ^-)z(x)z(x' )dxdx' }<?z . (2. l6) 

A 1 X ** X X 

Here the boundary condition on z(x) was taken into account: z(x) -► 0 for 

|p| oo and on the boundary of the configuration space. The second term 

in the curly bracket vanishes because 9<j>/3q = - . (Tnis illustrates 

that the sum of internal forces vanishes due to the action-reaction law.) 

The first term also vanishes if there is no external force (9H^/8q = 0). 

Otherwise it expresses the total external force. Therefore we can see in 

(2. 16) the generalized conservation law of momentum. It is needless to 

say that total momentum f npF dx is infinite in the limit N = nV -+ <*> but 

X 1 

it can be meaningful if we consider the average momentum per particle or oer 
unit volume. The other conservation laws are quite similarly examined, 
except for the energy c ise in which care should be taken of the 
fact that the effect of the term with a second order functional-derivative 
in (2. 9) appears for the part of the interaction energy including the ex- 
pression (2. 15 )* which is a quadratic functional in z. 

Finally we add a ccnment on the other functional formalism equivalent 

b 

to the BBGKY hierarchy by Nakayama and Dawson. In this formalism the 
characteristic functional plays a main role of describing the state accord- 
ing to the theory of Hopf 10 for turbulence mechanics. There the evolution 
of an ensemble is considered, all members of which are microscopic states 
evolving according to the Klimentovich equation 0 starting with various 
initial positions and momenta of particles given in the one-body phase 
space. The basic equation, i.e. the Hopf equation for this ensemble is 
somewhat different from (2. 2); in particular it is remarkable that it has 


0 $ . : >f • =’V/' ;s : *farfr:y ]**?%*; 
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no term with a coefficient quadratic in y(x). Correspondingly the cha- 
racteristic functional itself is different from the state functional 
(2. l). A difficulty with this formalism, however, stems from the fact 
that every solution of the Klimontovich equation does not represent a 
microscopic state. Hence all solutions of the Hopf equation are not 
realistic, even if they satisfy the general conditions on the character- 
istic functional. For the solution to be realistic, the members of the 

ensemble should be restricted to Klimontovich ’ s special class of singular 

8 

functions in terms of superposed delta functions. This limits the class 
of allowed characteristic functionals severely, and it would not be easy 
to formulate this restriction in a simple manner. Another disadvantage 
may be seen in the difficulty in finding an appropriate simple-form initial 
characteristic functional which can be used in practice; for the state 
functional it "s easy, as described in ref. 1 and will be shown in the 
next section. 


3. Master Equation Approach 

The basic equation (2. 9) is very close in form to the Fokker-Planck 
equation which describes a Markov stochastic process. The first term of 
the right-hand side is the analogue of the so-called friction term and the 
second resembles the diffusion term. However, there is a difference in 
principle between (2. 9 ) and the Fokker-Planck equation; the former is 
time-reversible as was mentioned already, but the latter is time-irreversi- 
ble. This reflects the fact that the coefficient function [ ; z(x)z(x T )] 
inside the second-order derivative is not positive-definite. This 
similarity and difference may contain a key to solving the noted 


historical question: how to bridge the microscopic rerersibility and the 




10 


macroscopic irreversibility in the system evolution in $ general way. 

Now if the coefficient function is considered as the sum of the 
positive-definite and negative-definite parts at each local point in A^, 
the former plays a role of creating an irreversible process, while the 
latter that of destroying it or creating an anti -irreversible process so that 
both roles may offset each other to cause the reversible process exactly. 

So, in order to extract a purely irreversible process from the dynamical 
process, evidently it is necessary and sufficient to retain only the 
former part of the coefficient function, so that (2. 9) reduces to a real 
Fokker-Planck equation. This would be the simplest universal way of 
introducing the positive time-arrow in the system evolution. There is no 
further assumption in this procedure, such as weak interaction, diluteness, 
etc. Therefore its application will be free from any such restriction on 
physical conditions. Let the Fokker-Planck equation thus obtained be writ- 
ten as 


ifi. = 

at 


-• f x 


I— {Qz(x)p}dx + 


<$z(x) 


If f L 

2 X J X 6 z(x) 6 z(x’ ) 


•{p[4>; z(x)z(x* ) ]p jdxdx' 
(3. 1) 


where the symbol p has been replaced by p to distinguish its approximate 
nature due to the above procedure and P is the operator on the coefficient 
function to make it positive-definite in the way just described. There can 
still be many ways of constructing a Markov process starting from (2. 9); 
for example, P[ <J> ; z(x)z(x f )l plus any nonnegative function can also be 
positive-definite. But it would be natural that (3* l) is the most appro- 
ximate to (2. 9) among them, since in this modification the first term is 
invariant and the second term was mi nimumly changed to make a real diffusion 
term. Then, the difference |p - pj would be minimum as a whole. Thus 
(3. 1) may be called the master equation in the sens.j that it masters 
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the irreversible system evolution which is the most loyal to the basic 
dynamical rule. 

By the property of the Fokker-Planck equation it is possible to inter- 
pret p as the probability density in the space A^. The condition (2. 11) 
plays an important role for this interpreration; that is 

/ p<5z = 1 . (3. 2) 

A 1 

If the diffusion term is neglected in (3. l), it can be found that the 
equation governs the time-evolution of the probability on an ensemble of 
the trajectories in the space A-^ which develop from various initial values 
according to the Vlasov equation: 3z/3t = Qz. This fact is obvious from 
the perfect analogy with the Hopf equation in turbulence mechanics * 0 if 
we go back to the equation for the characteristic functional ^ which 
corresponds to p (defined in the same way as (2. 7 )). In this case, we 
really deal with nothing but a so-called turbulent field which is basically 
governed by the Vlasov equation instead of the Navier-Stok.es equation. 

Here we find a complete coincidence with the other statistical-dynamical 

6,7 

research based on the Vlasov equation mentioned in Introduction. If 
the diffusion term in (3. l) is included, the coincidence breaks down and 
the motion (expressed by z) is affected by some internal random force im- 
plied by p[<fr; z(x)z(x* )], as is well expected from the Langevin equation 
for Brownian motion; so that wnat we may call a turbulence will be more 
irregular and more random. For this case there is the analogy with the 
Novikov equation in the turbulence mechanics with random force action 
(though nere is a slight difference in situation in that the random force 
in the Novikov equation is not internal but externally given independently 
of the field z). From all these facts, it is quite reasonable to in- 
terpret z(x) as the stochastic particle-number density in X normalized 
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in the sense of (2.13). Obviously z(x) is not a macroscopic observable 
(as is also known from the relations (2. Ik) and (2. 15), nor is it a 

g 

microscopic density such as considered by Klimontovich since it can be a 
regular function; but may be understood as a fictitious, mathematical 
working field with an intermediate property. 

Let us examine if (3* l) is consistent with the conservation laws in 
spite of its modification by P. It is known from the discussion of the 
previous section that the term with the second-order functional-derivative 
in the basic equation has no effect on the conservation laws escept for 
conservation of energy. Accordingly it is evident that the conservation 
laws of mass, momentum and angular momentum hold also for (3* l) . For the 
energy case, however, we have 

2 

-U( f nH (x)F (x)dx + f f B- <t >( q - q.' )F (x, x' )dxdx' ) 

3t J x 1 1 'XX 2 2 

2 

= — f f f z(x)z(x’ ) ]dxdx ! p6z , (3. 3) 

2 ' A i'x'x 

/ 2 2 . 

which can not vanish in general* The right -han^. side is 0(n X ) 

(X: representative order of magnitude of 0), but the error will accumulate 

with time. Therefore it is necessary to correct this error in our master 
equation some way. This can be made by confining the space for z into a 
manifold in which the conservation of energy is strictly insured. Indeed, 
since the total energy found on the left-hand side of (3- 3) is rewritten 

in terms of (2. 14) and (2. 15) as 

r % 2 

*' A p U nH i z (>0ix: + / / <fz(x)z(x f )dxdx f }6z , 

i A A 
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we can impose the condi cion on z: 

2 

f nH n z(x)dx + / / — 4>z(x)z(x’ )dxdx r = const. (3. *0 

'X 1 XX 2 

to guarantee the constancy of the total energy. Equation (3. *0 pre- 
scribes a subset of A 1 , which may be called a quadratic hypersurface. 

It is thus essential that our basic equation (3. l) should be re- 
defined in such a Riemannian space as prescribed by (3. 4). There is no 

special difficulty in dealing with the Fokker-Planck equation in a 
12 

Riemannian space. As is well known, a Fokker-Planck equation governs 

a generalized Brownian motion and hence we may imagine from our equation 

a generalized Brownian motion or random walk in the Riemannian space. 

It is known that there is a unique steady state of p, and p tends 

13 

asymptotically to this state irrespectively of any initial condition. 

It is unfortunate that we cannot at present have any explicit form of 
the steady-state solution to be compared with the equilibrium canonical 
distribution for the Gibbs ensemble. However, some similarity in 
situation may be expected. Namely, the final asymptotic state is the state 
in which all Brownian particles have fully diffused over the whole 
Riemannian space = the quadratic hypersurface (3. h) that corresponds to 
the constant-energy shell in the grand (6N-dimensional) phase space; and 
then every point in that space is realizable with some measure irrespec- 
tive of an initial state, just like a grand phase point in the energy 
shell in the equilibrium state. We note that in our theory, the concept 
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of entropy is not necessary to conclude the one-directional evolution of 
the system, since we have already the uniquely given asymptotic state. 

Or 

The initial condition imposed on p is easily formulated, if the 
initial correlations among particles are neglected. Namely, on starting 
from the initial state functional 1 

i> - exp{i/^y(x)F c (x)dx} (3. 5) 

A 

where F Q (x) is the given initial one-body distribution function, we have 

p = 6 [ z(x) - F Q (x) ] (3-6) 

by means of the relation (2. 7 ). (See ref. 1 for the delta functional 
6 [ ].) Naturally, here ve rewrites p as c . This means that our random 
walk begins from the single point in the Riemannian space which is given 
as F q (x). Then, the total of the Brownian trajectories beginning from this 

point provide all the physical information of the irreversible system 

f \ J 

evolution. Indeed, as is known from comparing (2. l) and (2. 7), all F g 
(the symbol F g means the approximation to F g in the sense of our master 
equation approach) are given as the s-th order correlation function of the 
stochastic field z(x); 

F s (x x , x g , t) = / z(x 1 )* *z(x s )p(z, t) 6 z (3. 7) 

where the integration is taken over the Riemannian space in and all 
physical quantities are related to F . It may be further noted that since 

the fluctuation of a macroscopic quantity is closely related to the behavior 

2 — 
of F£ , an observable macroscopic turbulence should be related through 

(in ( 3 . 7 )) with the turbulence in z which is strong enough to have a correlation 
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length over a macroscopic scale. This is an interesting situation which opens 
a new way of statistically pursuing a turbulence in plasma. 

In order to compute p in (3- 7) f° r a given M, it is necessary to 
solve the Fokker-Planck equation with many variables. A practically useful 
formula ior solving the Fokker-Planck equation is presented in the next 
section for the simple case where the Riemannian space can be approximated 
by a Euclidean space. This formula tells us how to follow the generalized 
random walk with the aid of a high-speed computer to constitute the 
solution of our initial -value problem. It will turn out that the principle 
of the calculation is nothing but the Monte Carlo quadrature. This method 
of solution is in principle available also to the (not simple but) 
general case, for which we should consider the random walk on the Riemannian 
space. Indeed, since such a Riemannian space is always approximated by 
the tangential hyperplane locally, the infinitesimal random-walk motion 
around the tangential point can be followed by quite the same way as given 
for the simple case, and continuation of this process would complete a chain 
of random walk on the hypersurface. 

Finally, it is interesting to point out that a new hierarchy similar 
to the BBGKY hierarchy can be obtained from our master equation. Let us 
multiply (3* l) by z(x) and integrate it with respect to 6z over the 
Riemannian space, taking into account (2. 10), (2. 14) and (2. 15). Then, 
we have 



(3. 8) 
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in the limit M -* 00 . Here use was made of partial integration with respect 

to <5z and the closure property of {s (x)}. This is formally the same as 

J 

the first equation of the BBGKY hierarchy. However, if we multiply (3. l) 

by z(x)z(x*) and integrate it to obtain the second equation which governs 

Fg, we find easily that the equation obtained is not the same as the 

second equation of the BBGKY hierarchy because now the second term of 

(3. 1) comes into play to give an effect of the modification due to the 

operator P. The same is true for the higher-order equations of the 

hierarchy. Thus, Fg, F , ***cannot be time-reversible, so that F^ is 

also time-irreversible. This is an essential difference between (3. 8) 

and the corresponding equation for From the viewpoint of the new 

hierarchy, it is no puzzle that the proper irreversible approximation to F2 

in (3. 8) as seen in Born and Green* work gives rise to the Boltzmann 

15 ~ 

equation with irreversibility. One may see a coarse-grained aspect of F^ 
from this fact, but we will come back to this matter in the final sub- 
section of Sec. 4-. 

4. Application to the Weak Interaction Case 

By the weak interaction case we mean the case where the total interac- 
tion energy between particles is far less than the whole energy of the 
system, such as the case of a gas or sometimes a plasma. In this case, 
the Riemannian space reduces to the fixed hyperplane in which 

satisfies 


n / H (x)z(x)dx/V = e 
X 1 


(e : const . ) . 


<fc. 1) 
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*r- 




This constitutes a subspace of A^, which we may call B. Or problem thus 
reduces to solving in general the master equation (3. l) in B. Here we 
give the explicit formulation of the operator P and next describe the 
practical program of how to calculate a physical quantity on the basis of 
our functional random-walk model with the use of a high-speed computer in 
mind. 

(i) Euclidean coordinates in B 

To represent a point in B in terms of a set of the Euclidean coordi- 
nates, we derive an orthonormal function set in B, starting from {s (x)}. 

J 

The limiting process V -► » is put off until we arrive at the final formula. 

Now in the representation: z(x) = £bjSj(x), (2. 15 ) defining 
is given as 


l b jJ x s j ( x)dx/V = 1 


A 


(h. 2) 


To find a rotation of 


^ su 


such that one of the basic vectors becomes a normal 


to the hyperplane A , the transformation matrix (t. .) such that 

1 ij 


- 1 


(u. 


kj 


or 


M) = 


1 1 


(U. 3) 


is introduaed. {b£*b are the new coordinates and the normal coordinate to 
is taken as bjj^b If we put b^ 1 ) = 0 except for k = 1 and b^ 1 ^ - 1, 
then (U. 3) gives 


* _ k X 

*U 


[ !(/ s. (x)dx/V} 2 ] l/2 


/ 8, (x)dx/V 
Y J 


» 


(b. k) 
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since {b } in this case are nothing but the direction cosines of a normal 
J 

to A . This fixes a part of elements of the matrix (t ). The other 
1 i j 

elements, however, can be arbitrarily given except under the condition 
that (t ) should be an orthogonal matrix; this arbitrariness corresponds 
to the freedom of space rotation within A . 


In the new coordinate system the hyperplane can be expressed by 

b^ 1 ^ (a certain const.) Hits constant, which is the distance 

M 

between the hyperplane and the origin of R , is found by the substitution 

of ( k , 3) into (1+. 2) together with b^ ^ * 0 except for k « 1. Namely 

k 

0) 


v - 1/{ \ y x s j (x ' dx/v} 

or l/[£{/ s (x)dx/V} ] ^ . 
j X 0 


(U. 5) 


Thus, a function in A^ can be expressed as 


(x) = l bl 1 M 1 )(x) + c (x) 
j>2 ° J 1 


(l*. 6) 


where 


sj^(x) = l t s (x) , 
J jk k 

c (x) = b^ 1 1 ^(x) . 

1 lo 


(fc. 7) 
(U. 8' 


In quite a similar way, we can find a rotation of such that one 

of the basic vectors becomes a normal to the hyperplane B. This time we 

. ( 1 ), 

note that the transformation matrix (t. ) is of M - 1 dimensions with 

ij 

the suffix conveniently begining with 2. Corresponding to (k. 1), we have 


fH s^(x)ix/V 

t' = '—LA j>2 

[ l (J H s^ 1) (x)dx/V} 2 ] 1/2 

k>2 * 1 k 


(k. 9) 
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on account of (4. l). The distance between B and the origin of is 


given as 


e - n/H^c 1 (x)dx/V 


20 ■>( I ILh , s< I) (x)dx/V} 2 ] 1/2 

j>2 X J 

As a result , a function in B is expressed as 

z(x) * l b^s^(x) + c (x) + c. (x ) 

,1>3 J J 


(4. 10) 


(4. 11) 


where 


'(*) = I t j^ s k 1 ^ x ^ » 


(4. 12) 


c 2 (x) = . 


(4. 13) 


Hence we know an orthoncinal function set in B is given by (U. 12) 


s| 2 )(x) = l I t k iS i (x) , J > 3 

1 k >2 Jk 1>2 kl 1 


(4. 14) 


( 2 ) 

and {bj ) is a set of the Euclidean coordinates based on it; but note that 

the space spanned by (4. l4) is not equal by itself but parallel to 

B, as is seen from (4. 11 ). 

(ii) The explicit form of p 

M —“P 

In the frame of f( now obtained, z(x) in the master equation (3. l) 
should be replaced by 


z(x) * £(x) + c (x) + c„(x) ; 

X d 


(4. 15) 


here and hereafter a function with a bar on top indicates that it belongs 
to § = (spanned by (4. 14)). Differentiation is invariant to this 

replacement; 6/6z(x) * i/6z(x). In this frame, let -os consider the func- 


- 20 - 


tional quadratic form with the coefficient kernel function [<}>; z(x)z(x’ )]. 
After the orthogonal transformation in B to make it diagonal, it can be 
written as 


^ J^y(x)y(x f )[$; z(x)z(x* )]dxdx' = £ D^(z) 


(U. 16) 


where a* are the Euclidean coordinates of y(x); the (real) eigenvalue s D. depend 

j J 

on <j> and z. Thus, the minimum modification of (U. l6) to make it positive- 
definite is simpDy to neglect D (z) when D (z) < 0. This leads to the 

j J 

formula : 


Vx* (l * (x,)PU; z (x) z (x' )]dxdx' = l a^ 2 D^(z) (4. 17) 

where 

D (z) = D (z) , P > 0 ; 

J J J 

= 0 , D i 0 . (4. 18 ) 

J 

Although this means just ,f completely prohibiting the anti-irreversi- 
ble process to happen”, as already described, an interesting relation with 

the coarse-graining idea will be explained in the final subsection. 

, ^ 

(iii) General solution for p 

% 

A general solution for p can formally be constructed in the form of 
a repeated multiple integral by the propagation kernel methodl’ 16 That 
is, 

p(z, t) = lim / ••/ P (z L /z L-1 )* *P (z 2 /z 1 )p(z 1 , 0) n 6z k , (h. 19) 
At+o B B At k=l 

where z » z L , At = t/L, and the superscripts indicate the order of the 


^ i < >ll< ! P ^ • 1 i >ff^^ iw i| ^ ^n^iiyo^ >» ■ - *•.. . . . . .*rq w »< , «** 
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time-subintervals. The infinitesimal kernel P is explicitly obtained 

At 

from ( 3 . 1 ) as 

P (z k+ 1 /z k ) = / exp[ i / y k (x){z k (x) - z k+ 1 (x) + AtQz k (x)}dx 
At g X 

- ^ / / y k (*)y k (x')PU; z k (x)z k (x')]dxdx']6y k , (1*. 20) 

2 XX 

and p(z\ 0) is the initial condition on p. It should be noted that the 

% _k Jr / \ \ 

arguments in p and P . are actually z rather than z through (4. 15) in 

At 

the present frame of B. The present notation of arguments is entirely for 
convenience. 

The expression (4. 20) may be rewritten in terms of the Euclidean 
_k, x 

coordinates for y (x) ; 


y k (x) = l a k *s k *(x) , 


(4. 21 ) 


where {s k *(x); J i 3) is the orthonormal function set in B which is related 

J 

with (4. 14) by the orthogonal transformation described in the preceding 
subsection. By noting (4. 17) > we have 


P At (zk+1/zk) = l ak * {6k * " bk+1 * + At/ v S k *(x)Qz k (x)dx} 


>3 J J i 


Ji3 


’X j 


- A| l ( a k *) 2 J (z k )] n {de k */( 2 ir ) 1/2 } , (4. 22 ) 

2 J>3 J 


J J>3 J 


where 


k* . v# \r 

6 = / S ( x ) z (x)dx 

J X J 


(4. 23) 


and 



b k+1 = / s k *(x)z k+1 (x)dx , 
J & J 


which is further calculated as 


(4. 24) 
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p At u k+ V) 


J Utfyz*)) 1 /* 


{0 k %At/ e k *(x)Qz k (x)dx-b k 

exp[ J fr .. J A 


k+l* } 2 


2AtBj(z k ) 


-J • 


(^ 25) 


Here we note that the j-th factor with Dj = 0, if av. reduces to the delta 

function 6(e k# + At/ s k *(x)Qz k (x)dx - b k+1 * ) . 

j X o « 

p in (4-. 25) is a generalized Gaussian measure factor with respect to 
At 

wt* 

{b. }, so that it may be interpreted as the probability density for the 

J 

k+ i* k 

stochastic quantity {bj } . Since this probability depends on z and 
since {b k+ ^ } represents z^ +1 in the manner of (4. 24), it may be understood 
as the transition probability that a random-walker in the space B (or B) 


moves from the point z k (or {b^*l ) to the point z k+1 (or {b k+1 }) in a small 
time At, A great difference in our random walk from the simple Brownian motion 

k 

is seen in the complicated dependence of P^. on z , i.e. the starting point 
at each instant. As a consequence, the measure p^z established by (4. 19) 
together with (4. 25) is far more complicated in shape than the Wiener 
measure. However, it is worth noting that the present master equation 
approach has the great mathematical advantage over the original state func- 
tional approach^ - in that the existence of the integral t = J*_ exp (ijyzdx) p'AT 

~ ^ 17 

is guaranteed by establishment of the measure p bT even in the limit M -* ® 1 . 

Thus, it is concluded that all Fs converge as M - « . 

(iv) The Monte Carlo quadrature 

It is now obvious that if (3* 6) is adopted for p*(z^, 0), (4. 19) ex- 
presses the probability of finding the random walker, who started from F q (x) 
in B, at z(x) after a continuous chain of random-walk during the time t. 

The Whole integration in (4. 19), however, is analytically impossible. It is 
rather fortunate that the first integral with respect to z~ can be calculated 
as 




/. P A t (z2/zl) ^ (zl » 0)6zl = /- p At (z2/zl)6[5l(x) + C 1 ^ X ) + c 0 (x) ‘ F 0 (x)]6z ] 


- P At (z /F o> * 


(4. 26) 


o 

But the next successive integrals with respect to z , z **, i.e. in the 

Euclidean coordinates {b'T }, {b^ }, ***'are involved enough to recommend 

J 0 

the Monte Carlo quadrature with the use of a high-speed computer; although 

then we should be satisfied with an approximate value to the true integral. 

It is nice for an importance sampling to be applied that the integrand 

itself constitutes a product of conditional probability densities. Then, 

2 * 

the first task to estimate the second integral with respect to [b. } is Just 

J 

r 2 *t 

to sample values for [b. j out of the ensemble with the probability dis- 

3 

2 2 3 2 

tribution P^(z /F q ) and bo insert them through z (x) into P^Cz'yz ). 

({s. (x)} is known as a set of eigenvectors w_oh the eigenvalues {d.(F )}. 

J J ^ 

—2 ^ 2 
Hence z (x) is solved for by (4. 24). Together with {D,(z )}, a new set of 

J 

2*/ r 2*^ 

eigenvectors {s. (x)J is found and hence [j3 . } is solved for.) Thus, the 
J J 

3* 

probability for [b. } is decided and so the same sampling process can be 
J 

done for f } , and further for {b } , . . . until we arrive at (b^ ^ , when 

J J J 

one round of the importance sampling for the whole integral (4. 19 ) 
finishes. The next task is, of course, the average over many rounds of 
sampling . 

If we wish to estimate F by the formula (3. 7), another sampling is 

s 

r L*, 

necessary for the integration with respect to Ib^ } . In this case we only 
have as the estimator a function z(x ) .. z(x ) made of the sampled values 

JL S 

r L*t 

for ibj J, so that the average over many such sample functions may give 

F (x.,, .., x , t). This is the principle of the Monte Carlo quadratuie 
S «L s 
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for calculating F . Since (4. 25) is a Gaussian probability distribution 
s 

for {b^ + '*' }, we can sample the values for {b^" 1 "^ } out of normal random 
numbers with proper variances and averages. It is interesting to note that, 
each round of the importance sampling simulates nothing but each chain of 
random-walk during the time t in a segmented way. In the limit when At - 0 
or L - ® the simulation becomes perfect. Also, it is expected that the 
other approximations included in the present procedure are improved with 
the values of M and V increasing independently. 

(v) The relation with the coarse -graining idea 

We may point out a slight relation of our master equation approach with 
the coarse -graining idea. The prescription of (4. 17) - (4. 18) may also be 
understood to limit the space for y(x) into a smaller subsuace than R M ~ 2 by 
neglecting all the j-th eigenvectors of the space when D. ^ 0. This means 
that y(x) has been coarse-grained by losing many orders of orthonormal func- 
tions which are not desirable for the irreversible description of the system. 

However, this coarse -graining procedure is not so fixed an operation as the 

15 

conventionally understood one , but is very flexible because it depends on z. 

If we proceed so thoroughly with the physically motivated coarse -graining 

idea that we artificially change the domain of integration in (4. 22) from 

R to the subspace described above, the transition probability in (^. 25) 

will be the product of only the j-th factors with D^(z k ) >0. Then, b^ +1 

with the same j as that of the neglected eigenvector in y (x) does not appear 

in P so that z k+1 (x) and hence z k+1 (x) would be coarse-grained in the 

k -k 

ssme way (depending on z ) as y (x). Hence we can coarse -grain the distri- 
bution function F (in a flexible sense) through the relation (3. 7). The 
s 

above discussion based on the artificial change of the domain of integration 
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is interesting in suggesting the physical meaning of our master equation 
approach to some degree, hut obviously such an artifice is strictly inexact. 
Therefore, we can only see by this discussion a roughly sketched, probable 
physical picture of our theory. All what we may say is that some very 
flexible kind of averaging process must have been introduced as a result of 
the operation P. 

5 . Conclusion 

A minimum modiiication of the functional formalism of classical non- 
equilibrium statistical mechanics was performed (on the basis of the cylinder 
functional approach) such that the basic equation turns into the master 
equation, which includes the irreversible system evolution. The general 
theory was applied to the case with a gas or a weakly interacting plasma. 

The method of solution proposed in Sec. 4 has no great technical difficulty 
but for the use of a high-speed computer. All calculations with the Monte 
Carlo quadrature will be rather simple for the machine. It may be said that 
the possibility suggested in ref. 1 of approaching non-equilibrium statis- 
tical mechanics by the Monte Carlo quadrature has been clarified in this 
paper. It is expected that this method will clarify many -unsolved problems 
in relation with the irreversible process of the many-particle system, even 
for the case with spatial inhomogeneity. If we axe interested in a steady 
state, it can be studied by examining an asymptotic behavior of the solution 
as t -* <* for an arbitrary type of initial condition. Together with the 
tangential approximation to the Riemannian space, the method of solution can 
extensively be used for the general strong interaction case with an additional 
care. In these respects the practical value of the present theory depends 
to a large extent on future works employing a computer. It is interesting to 
note that a similar type of functional integral (with the Gaussian measure). 
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which corresponds to a solution for the .iiaracteristic functional for the 

Burgers model turbulence, was recently calculated by the Monte Carlo 

18 

quadrature with promising success 
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